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Abstract: Rabies is one of the neglected tropical diseases that has persisted for centuries in Ethiopia, and it is 
endemic within and around Addis Ababa. For the purpose of studying the dynamics of the disease we propose 
a deterministic mathematical model with human, dog and livestock populations and formulated as a system of 


ordinary differential equations. Basic reproduction number R, and effective reproduction number R, are 
computed using next generation operator. Sensitivity analysis of R, shows the natural death rate of dogs ug, 
the annual birth rate of dogs 94, dog-to-dog transmission rate Bg, and disease induced death rate og are found 
to be the most sensitive parameters of R, . According to numerical simulations of our system rabies 


transmission will increase within and around Addis Ababa, and may peak in 2024 and 2026 in human and 
livestock populations respectively. Our simulation shows that 2596 vaccination coverage in livestock 
populations will reduce the future infection by half. This study suggests that a combination of interventions 
consisting of 6096 of vaccination coverage in dog populations, 1596 culling of stray dogs, and reducing annual 
crop of newborn puppies by 2596 will reduce the number of human and livestock infections by 7096, and the 
disease will be eradicated from the community. 


1. Introduction 

Rabies is a zoonotic viral disease that causes an acute inflammation of brain in human and 
other warm blooded animals. It is transmitted by saliva of an infected animal via bites or 
scratches. Once the virus enters the body through a skin opening, it travels via nerve tissues 
to the brain where the virus duplicates itself. During this time the host experiences a range 
of symptoms from fever to hallucinations, paralysis, and eventually death [2, 9]. Salivary 
glands are attacked after the virus duplicates in the brain, then the saliva becomes the main 
instrument for infecting other animals. 

Rabies is the most fatal of all zoonotic infectious diseases. Unfortunately awareness about 
the disease is low and itis often recognized as a public health issue in developing continents. 
Despite the interventions and scientific breakthroughs, rabies continues to be a dreadful 
communicable infectious disease [1]. More than 55 000 people die due to rabies per annum 
in more than 150 countries and territories. In low-income countries the rabies surveillance 
systems are often poor and underreporting occurs because of infected people dying at home 
[10]. More than 40% of the people who are bitten by infected (rabid) animals are children, 
which is explained by the higher tendency of children to play with animals [18]. 
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Rabies is transmitted among animals and to human through bites or scratches of a rabid 
animal. Many mammals can transmit the rabies virus, but in many parts of the world rabies 
is spread through infected domestic dogs. More than 3 billion people in developing countries 
are exposed to dog rabies [9]. Dogs are the source for almost all types of human and livestock 
infections [5]. The saliva of an infected dog capable of transmitting the virus through 
scratches or bite to another animals such as another dog, cat, mule or cattle which can 
become rabid as well. The infected animal will demonstrate all symptoms and eventually die. 
Besides its effects on humans directly, rabies also puts a significant burden on the livestock 
population in Africa and Asia. This can be directly reflected by its economic effect on rabies 
endemic areas. Africa and Asia lose US$12.3 million annually because of deaths of livestock 
due to rabies [8]. 

Rabies is highly endemic to Ethiopia. It was estimated that 10, 000 people die due to rabies 
annually, which makes Ethiopia to be one of the worst affected countries in the world. Low 
awareness and lack of accurate quantitative information on rabies both in human and animal 
populations prevent people from applying effective control measures. In most cases 
individuals who are bitten by rabid dogs go to traditional healers, which interfere post 
exposure prophylaxis [8]. 

In most cases rabies transmission is very high in urban places because of the high number of 
domestic dogs [9]. The major cause of spread of rabies in these regions is urbanization. It is 
estimated that the dog to human population ratio in Ethiopia is 1:6 in urban and 1:8 in rural 
areas. The number of dogs in Addis Ababa is estimated to be between 150 000 to 200 000 
[14]. Despite vaccination and other control measures, rabies has persisted for centuries and 
it is reported from all regions of Ethiopia [15]. The first occurrence of rabies in Addis Ababa 
was recorded in August, 1903, and it was known by its traditional name called mad dog 
disease [13]. Addis Ababa and its surroundings are the endemic parts of the country. It was 
reported that around 2000 people are bitten by dogs annually within and around Addis 
Ababa [3]. A retrospective record review from 2001 to 2009 show that 386 human fatal 
rabies cases occurred, with an annual range of 35 to 58; From this, 42.72% were children 
under the age of 14 [5]. 

Mathematics has played a very important role in understanding and controlling the spread 
of infectious diseases, and it is a powerful tool for analyzing and predicting the dynamics of 
phenomena. It also helps medical professionals to organize their thinking [6, 11]. 

Some studies have been done using mathematical models to analyze the dynamics of rabies 
transmission among dogs and from dogs to humans, but none of them incorporate livestock 
populations which are the livelihood of poor communities. 

An SEIR (Susceptible-Exposed-Infected-Removed) standard model was developed by Addo 
[2] to determine and predict the spread of rabies among dogs in Bongo District Ghana. Both 
SEIR model with vaccination and SEIR model without vaccination were formulated with 
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ordinary differential equations. The reproductive ratio without vaccination was determined 
to be greater than one which shows the virus will be endemic, and less than one with 
vaccination, which shows that the disease dies out. The study also applied sensitivity analysis 
to the model by using different numbers of infectious dogs and vaccinated dogs. The study 
determines the reproductive number, R, of rabies transmission decrease as vaccination 


introduce to the model. In addition, the model shows that rabies transmission can be 
decreased by the strategy of keeping dogs confined within their household. 

The other deterministic SEIR model was developed by Zhang et al. [19] to analyze and 
control the transmission of rabies among dogs and from dogs to human in China. Both dogs 
and humans were included and classified into susceptible, exposed, infectious, and 
recovered classes. They first simulated human rabies from 1996 to 2010 using the data 
reported by Chinese Ministry of Health, and the numerical simulation they got significantly 
supported the data. They also estimate the basic reproductive rate R ~2 for rabies 


transmission in China. Sensitivity analysis of R, was performed in terms of the model 


parameters and compare the effects of culling and immunizing of dogs. Their results show 
that reducing dog birth rates and increasing dog immunization coverage rates are the most 
effective methods in controlling human rabies infection in China. They recommend that 
culling of dogs can be replaced by immunization of dogs. The reason is that in the process of 
culling of dogs human community can be disturbed. Additionally culling can cause increased 
movement of infected dogs to less infected areas. Their model predicted that rabies 
transmission in China will decrease for the coming 7 to 8 years and it will peak again in 2030. 
SEIV (Susceptible-Exposed-Infectious-Vaccinated) model was formulated by Hou et al. [7] 
for the transmission of rabies among dogs, and from dogs to human in the context of 
Guangdong province of China. In their model domestic and stray dogs were taken as different 
groups, and the model was governed by twelve differential equations. In their study the 
sensitivity analysis ofthe reproduction rate was determined in terms of various parameters. 
The reproduction rate was R, ~1.65 which is less compared to [19]. According to their 


results the recruitment rate of domestic dogs, the number of stray dogs and the valid time of 
immunity play a very important role for the transmission of rabies. 

In this study we propose a mathematical model for the transmission dynamics of rabies from 
dogs to both humans and livestock in the context of Addis Ababa and its surrounding areas. 
The model is based on SEIR type and domestic dogs infect both human and livestock 
populations. We have used the improved model of Zhang et al. [19] to incorporate livestock 
populations. Though livestock populations are highly affected by rabies virus, they have not 
been considered in the above studies. 


4/23 
EGA, LUBOOBI, KUZNETSOV AND KIDANE 


2. Materials and Methods 
2.1 Model Formulation 

In this chapter we formulate SEIR (Susceptible-Exposed-Infected-Recovered) model of 
rabies for human, dog and livestock populations. We categorize the human, dog and livestock 
populations into susceptible, exposed, infected and recovered groups. Susceptible groups 
have no disease, but they are likely to be infected in case of contact with rabid dogs. Exposed 
individuals are those who contracted the virus via bites or scratches, but still they have not 
shown symptoms. Infected individuals are those who develop clinical symptoms and they 
are unlikely to recover due to the nature of rabies. The recovered classes are those who 
recovered through vaccination before they reach infectious stage, whereas the rest get 
infected and die eventually. 
All the parameters of the model are positive and they are introduced in Table 2.1 
Table 1: Description of parameters. 








Parameter Description 
9939 The annual birth of human, dog and livestock populations respectively 
E E per annum 
— Death rate due to rabies for human, dog and livestock populations 
PET respectively 
Donc The loss rate of vaccination immunity for human, dog and livestock 
EOD populations respectively 
Hs Las Hg Natural death rate of human, dog and livestock populations respectively 


The rate of infectious dogs infects susceptible human, dog and livestock 


Po Bi B, 


populations respectively 


Py PaP, The incubation period in human, dog and livestock populations 


resnectivelv 
Vaccination rate of exposed human, dog and livestock populations 


Pr» 0, , P 


respectively. 


The human population is grouped into susceptible, S,, exposed, En, infectious, I}, and 
recovered, Rh. Individuals are recruited to susceptible class by birth at 2n to susceptible 
human. A susceptible man bitten by a rabid dog becomes exposed. If post-exposure 
treatment is not given the person become infectious and dies since there is no recovery at 
infectious stage. 

The dog population is divided in to susceptible, Sy, exposed, E4 , infectious, I4 and 


recovered, R4. Individuals are recruited into susceptible class by birth at 2,. For susceptible 
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dogs vaccination is applied at a rate of g,. This is because it is the dog population which 


infects both human and livestock populations. An exposed dog moves to the infectious class 
directly since the community cannot observe which dog is infected as many dogs are very 
mobile around the city. 

Livestock populations are also divided into susceptible, S;, exposed, E;, infectious, I; and 


recovered, R;. Individuals are recruited to susceptible class by birth at 9. Members of the 


susceptible class contract the disease from dog bites or scratches. Livestock which are bitten 
by a rabid dogs become exposed. If post exposure treatment is provided the individual moves 
to recovered class before reaching the infectious stage. 

Our model is developed based on the following assumptions. Susceptible populations are 
recruited by birth at 9; Rabies transmission among humans, among livestock and between 
human and livestock was ignored due to negligible cases; An individual who is bitten or 
scratched by rabid dog becomes exposed; Individuals in each group have equal natural death 
rate; Populations are homogeneous, that is each individual has equal probability of being 
bitten or scratched by a dog and thereby contracting the diseases; Once an individual reaches 
to infectious stage there is no recovery and death is 10096 certain. 

2.2 Model Compartment and Dynamics 
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Figure 1: Flow diagram for rabies transmission among dogs and to human and livestock in 
which the parameters are as defined in Table 1. 

Using the assumptions, definition of variables and parameters, the model flow diagram 
which depicts the dynamics ofrabies transmission among dogs and from dogs to both human 
and livestock is shown in Figure 1. 
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The parameters of the model are positive. 9. where i = h,d,l represents the annual birth of 


dog, human and livestock populations respectively. Exposed populations of human and 


livestock can recover through post exposure treatment. The parameters p, where i = ^,d,l 


represent the latency rates of human, dog and livestock population respectively so that = 
i 
where i=h,d,/ are the corresponding incubation periods. 
2.3 Model Equations 
Based on the assumptions and interrelation between the variables and parameters in Fig.1 
rabies transmission dynamics can be described by using ordinary differential equations. 


dS 
— = 3, + O,R, — B,I,S, — 1,8, 


























dE 
205 = B,1,S, zj Co, tu, 9,)E, 
d Human (1a) 
dl, 
= p,E, - (4, +o, 
dt 
dR, 
dt = =0,E h - (0, F Hy )R, 
dS, 
= Q, to,R,-(u, +P; + B,I,)S, 
dE, 
= BS ily - Co, + Ma) Ea 
Dog (1b) 
= p, E, - Ua t+ 04), 
"m = PS4 — (La + @,)Ry 
dS, 
=9+0,R,-61,8,- “5, 
dE, 
po = BI,S, - (p, + uj + QE; 
Livestock (1c) 
dl, 
— = pE, - (uj +0), 
dt 
dR, 
zi = QE, — (0, + M )R, 


The total human dog and livestock populations are N, (t), N (t) and N,(f) given by 
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N,(D=S,(D+E,(DO+L(O+R, O 
N,(t)=S,() +E, (0) +L (OD +R, (0) 
N(D=S,(D+E(O+L(O+RO 


Therefore adding each of the differential equations of system (1) of human, dog and 
livestock populations will give us: 


dN, 

—£ =Q —u,N,—o,l 

di h THING» nån 

dN 

—— = 9 — tg Nga = Fala (2) 
dt 

dN, 

— = 9 -yuN,-o0;l 

dr 1T HIN, 15i 


where N,, i 2 A,d,l is the total of human, dog and livestock populations at time t. 


Invariant Region 

The model system represented by (1) dealing with human, dog and livestock population will be 
analyzed in the feasible region, and all state variables and parameters are assumed to be positive 
for all t 2 0. The invariant region will be obtained through theorem 1. 

Theorem 1: 











All solutions of the system (1) are contained in the region  €(17 and D=®, U®, UD, 
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Proof. The model system (1) was grouped in to Human N, Dog N, and livestock N,. Such 
that 
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and © is the positive invariant for system (1) 
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Thus P=, UD, UD, el; xU xu). 























From this it is sufficient to consider model system (1) in the region 6, and it can be shown 
to be positively invariant The model can be considered as epidemiologically and 
mathematically well-posed. 
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3. Model Analysis 
3.1 Disease Free Equilibrium Points (DFE) 
To find the disease free equilibrium points we set the right hand side of equations of system 
(1) equal to zero. In the absence of attack or in the absence of rabies the following 
compartments will be zero. 


E,=R,=l,=Ej=L=E=R=1=0 





then the disease free equilibrium (DFE) & will be 
£ = ($5,0,0,0, $,,0,0, R, S,,0,0,0) 





viesta us By tg + Og) pk, = P19 and so 
Hn My (Ma + 9, + 0,) My (Ma + 9, 04) Hi 
For the dog population in the case of disease free equilibrium points Å, cannot be zero 


because susceptible dogs which are vaccinated transfer to recovered class. Therefore the 
disease free equilibrium points of system (1) is exist and given by: 
äi -( 2.0.0.0, eae Oi e 20 00.0] (3) 
H, Ha (My t Pq + Oz) Halla +P t0) 4 
3.2 The Basic Reproduction Number, R, 





The basic reproduction number œ, is a threshold parameter defined as the average number 
of secondary infection caused by an infectious individual by introducing in to a completely 
susceptible population. It is also called basic reproduction ratio or basic reproductive rate 
[6]. Basic reproduction number is very important parameter which helps to determine 
whether the disease spread in the population or it dies out. If R, <1, then on average an 
infected individual produces less than one new infected individual over the course of its 
infectious period, and the infection cannot grow. Conversely, if R, >1 then each infected 
individual produces, on average, more than one new infection, and the disease can invade 
the population. It is also crucial in the process of analyzing sensitive parameters which drive 
the dynamics of the disease and stability analysis of disease free and endemic equilibrium 
points. 

To compute the basic reproduction number it is important to identify new infections from 
all other changes in the population. We used next generation operator method proposed by 
[16]. We considered system (1) without vaccination i.e. o — 9 — 0. 


Let f(x) be the rate of appearance of new infection in compartment i, v; (x) be the rate of 


transfer of individuals out of compartment ; and v; (x) be the rate of transfer of individuals 


into compartment ; by all other means, and it is assumed that each function is continuously 
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differentiable at least twice in each variable. The disease transmission model of system (1) 
consists of nonnegative initial conditions together with the following system of equations: 


ž = F(x) = fl): i=1,...,6 (4) 
where v, =v; —V, 

We consider expressions in which the infection is in progression. These are 

EI ET, E and T. 


i = B,S,I, - Co, + 4 JE, 
dt 

dl, 

dt 

dE, 
dt 


dl 
—- p, E, - (4 + Ogg 


= p,E, (4, * 0,)H, (5) 








= B,S,I, (Pat uj) E, 


dt 
dE, 
— = p Sida- (P, + ME, 
dt 
dl, 
— = pE - (4 +07)1 
di PE, - (u * 01), 
By rearranging equations of system (1) without vaccination from infected to infectious class 
of human dog and livestock population with a system of equations given by (5) 


Let F bea non-negative nxn matrix and V be a non-singular N-matrix such that 


Xj Xi 


The point &, is the disease free equilibrium point in (3) without vaccination. 














where 
f BSL |v, | I Go, +h) En 
f; 0 Vy (Hy, +O, Ly - PrEn 
gu f i PaSala id im (Pat Mg) Eg 
1 ds 0 V4 (4, + 04), - p,E, 
fs BSI, Vs (o, + ME, 
fs 0 (6) | v | | Cu, tol, -pE | (7) 


By considering the classes in which infection is on progression and using the linearization 


technique. The Jacobean matrices of fand v at the disease free equilibrium point &, are 


given by: 
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o 00 h 9 0 
Hn 











Of (Ey) o 0 0 /* gp 
pale we m (8) 
"n 
0 00 0 0 0 
0 00 PÅ, 0 0 
H; 
0 00 0 O0 0 
[po + Hy 0 0 0 0 0 | 
— n Hy tO, 0 0 0 0 
y= iE) — 0 0 pat, 0 0 0 (9) 
Ox, 0 0 =p, fig >; 0 0 
0 0 0 0 Py Ui 0 
| 0 0 0 0 -D hte, | 





2 
Solve for v^ and multiplying it with F will give us a matrix ry-! = 2 (59) | | ev; 2l 
X 


Xj 


Jj 
Choosing the maximum eigenvalue in absolute terms that is the spectral radius ofthe matrix 


FV" gives us the basic reproduction number which is given by: 


= Pab 
Hy (Py + My Xu, +04) 





(10) 


In our assumption there is no secondary infection in human and livestock population, due 
to this all the parameters of R, in (10) are depending on dog population. This shows that 
targeting dog population in the process of combating rabies transmission is very important. 
Providing PEP for human or livestock population does not eradicate the disease from the 
community. 
3.3 The Effective Reproduction Number, R, 

The effective reproduction number is defined as the measure of average number of 
infections caused by a single infectious individual introduced in a community in which 
intervention strategies (in our case vaccination) is administered. The effective reproduction 
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number R, of system (1) is computed by applying the same procedure of R,. The spectral 


radius (dominant eigenvalue) of rv denoted by R, = p(FV ). 








R = D, PO, Gu, + 07) (11) 
usus + 9, t o, Xp, + Ha Ma to) 
Ó 
Since Ry = PaPsda , we can express R, in terms of R, as 
Hy (Py + uu Xu, +03) 
+@ 
_ Qu to) R (12) 


g (la * 9, + 09) 
The simulation for the basic reproduction number shows that rabies transmission is very 
high without any intervention. We simulated the effective reproduction number with 
different vaccination coverages and a combination of interventions consisting of 
vaccination, culling and controlling the annual crop of newborn puppies. The results show 
that as we increase vaccination ofthe dog population there is a possibility for the disease to 
die out. The simulations show that a combination of interventions are the best way to 
reduce Re less than unity. In Fig. 2, Rois without any vaccination, AK is the current 10% 
vaccination coverage, Ae is 60% vaccination coverage, Aes is 9096 vaccination coverage and 
Rea is combination of interventions consisting of 60% vaccination, 15% culling of stray dogs 


and reducing the annual crop of newborn puppies by 25%. 
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Variation of Reproduction Number with Exposure rate B å 
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Figure 2: Reproduction number for different vaccination coverages and combination of 
vaccination, culling and controlling newborn puppies. 
From Fig. 2 we observe that R,«R,«R,«R,«R, which shows that as we increase 


vaccination of dog populations the effective reproduction number decreases and becomes 
less than unity. It is very important to combine different interventions so as to facilitate the 
eradication of the disease form the community. 

3.4 Sensitivity Analysis of R, with Respect to the Model Parameters 


The aim of any mathematical epidemiology study is to understand the dynamics of a disease 
so as to control it by targeting some sensitive parameters. This can be achieved by 
performing sensitivity analysis based on the model parameters. The parameters of our 
model are given with their descriptions in Table 2. 

Sensitivity analysis tells us how each parameter of R, affects its result. This will help to 
identify which parameters are most sensitive for the spread of the rabies virus, so that 
appropriate control measure can be taken [4]. Sensitivity analysis can be done by computing 
the sensitivity indices of Re. With small variation of the parameters it is important to identify 
which parameters greatly affect Re. 
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Table 2: The parameter values of the model 
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Paramet Description Value (Year?) Source 
er 

9, Annual birth of humans 12 1980 Assumption 
0), Human loss of vaccination 1 Assumption 
P, Human incubation period 0.17 Zhang et al 
i Natural death rate of humans 0.016 Assumption 

fi Dog-to-human transmission rate 1.29 x10? Fitting 

Ø Vaccination rate of humans 0.54 EPHI 

1 

0, Disease related death rate of 1 Assumption 
9, Annual birth of dogs 2x10? Fitting 

0, Loss of vaccination immunity of — 0.5 Assumption 
p, Incubation period of dogs 0.17 Zhang et al 
li, Natural death rate of dogs 0.083 Assumption 
f, Dog-to-dog transmission rate 1.29x10> Fitting 

9, Vaccination rate of dogs 0.1 Assumption 
9: Disease related death rate of 1 Assumption 
9, Annual birth of livestock 2x10? Assumption 
5 Loss of vaccination immunity of 0 EPHI 

i livestock 
P; Livestock incubation period 0.17 Assumption 
Li Natural death rate of livestock 0.05 Assumption 
B, Dog-to-livestock transmission 1.18x10? Fitting 
9, Vaccination rate of livestock 0 EPHI 
Disease related death rate of 1 Assumption 

0, 


livestock 


The normalized forward sensitivity index is the ratio of relative change of a variable to the 
relative change in parameter. Ifthe variable is a differentiable function ofthe parameter then 


the sensitivity index is defined as follows: 
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Table 3: Sensitivity indices of R, 








Parame Description Value Sensitivity Indices 
i Natural death rate of dogs 0.083 -1.3781 
9, Annual birth of dogs 2x10? +1 
f, Dog-to-dog transmission rate 1.29x10° +1 
9, Disease related death rate of 1 -0.9259 
Pi Dog incubation period 0.17 0.3243 
9, Vaccination rate of susceptible 0.1 -0.1471 
dogs 
Oj Loss of vaccination immunity of 0.5 0.1268 
dogs 


Definition 1: The normalized forward sensitivity index of variable g that depends on 


parameter b is defined as: 
SUE 8 


8 — 
=>" 


(13) 


Since we have computed the effective reproductive number, R,, the normalized forward 


sensitivity with respect to the parameter b is given by: 





OR, b 
8 — ex 
Ys ob R, 
For example the sensitivity indices of R, with respect to Å, is given by: 
OR, 9 
8 =e yd =+] 
Ps 09, R, 


(14) 


By using the same notion, the sensitivity indices of the effective reproduction number given 


in expression 11 is computed with respect to all parameters embedded to R,. 


According to the sensitivity indices, the parameters of annual birth of dogs 9, and dog-to- 


dog transmission rate /, are the most positively sensitive parameters. This means 


increasing the parameters increases the effective reproduction number and vice versa. For 


instance increasing 9, by 10% will increase R, by 10%. Decreasing 9, by 10% will 


decrease R, by 10%. Dog loss of vaccination immunity 9, and dog incubation period 9, are 


less sensitive positive parameters. 


Natural death rate of dog u, and disease related death rate o, are the most negatively 


sensitive parameters. Increasing the parameters which have negative signs will decrease the 
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value of R, and vice versa. For instance Ve =-1.3781 means increasing the natural death rate 


of dogs by 15% will decrease R, by 18%. 


In summary, targeting the most positively and negatively sensitive parameters in the process 
of combating rabies will be most effective in reducing the transmission of the virus within 
and around Addis Ababa. 
4. Numerical Results and Discussion 

In this section we present the numerical simulations of the system (1). We used ode45 
MATLAB's standard solver for ordinary differential equations (ODEs). This function 
implements a Runge-Kutta method with a variable time step for efficient computation. Using 
the data reported from EPHI we first simulated human rabies fatal cases from 2008 to 2014. 
The data like the annual birth of dog populations, dog-to-dog transmission, dog-to-human 
transmission and dog-to-livestock transmission are not found from EPHI data, they are 
obtained by data fitting. 

In Fig. 3 the broken curve represents the data which was reported from EPHI and the smooth 
curve is the simulation of the system (1) of infected human populations. There is a good 
match between our model and the reported data. Our simulation further predicts that human 
rabies will increase for the upcoming thirty two years. The simulations are based on 
parameters which are presented in Table 2. The initial conditions of the variables are 
consisting ofthe data that are reported from EPHI, logical assumptions and data fitting. Thus 
S, (0) — 5x106, E,(0) =100, I, (0) = 38, R,(0) =2.5x 104, $,(0) 23x105, E,(0) = 8000, I, (0) 


—4000, R, (0) 25x 104, E,(0) = 2.5x 105, E,(0) —50, 1,(0) —5, R, (0) =0. 
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Figure 3: Comparison between reported data and simulation of system (2.1) for rabies 
infected humans in and around Addis Ababa from 2008 to 2014. 

In Fig. 4 the numerical result of the infected human population shows that the rabies virus 
will spread very rapidly in the coming years and it will peak in 2024. The bite of a rabid dog 
is the main reason for the transmission of rabies to humans. The increased number of dogs 
within and around Addis Ababa raises the number of stray dogs. Surveys indicate that in 
Addis Ababa, 1299 dog bites and/or scratches were reported in humans for the period 
September 2008 to August 2009. The majority of bites were made by stray dogs [12]. This 
indicates that there is a need for stray dog control strategy to reduce the number of human 
infections. 

From our sensitivity analysis, we found that annual dog birth is one of the sensitive 
parameters which controls the dynamics of the disease. A minor increase in newly born 
puppies increases human infection and vice versa. Fig. 5 shows a different figure of annual 
birth of dogs versus human infection. It can be noted that applying a strategy to control the 
annual birth of new born puppies is one of the most effective ways to reduce human rabies. 
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Figure 4: The dynamics of rabies in infected humans for the next 32 years in and around 
Addis Ababa. 
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Figure 5: The effect of annual birth of dogs 92, for human rabies infection. 


From the data fitting the annual birth of dogs is estimated to be 20 000. Replacing this value 
with a lower estimate of 18 000 annual dog birth significantly decreases human infection. 
This shows that targeting the annual birth of the dog population is very significant to reduce 
human infection. 
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Figure 6: Transmission of rabies in dog population in 40 years' time. 

The numerical solution for the dog population in Fig. 6 shows that based on the current 
condition of rabies transmission dynamics the disease will not perish. As the susceptible dog 
population decreases, exposed and infected populations increase. Dog-to-dog transmission 
is higher than dog-to-human or dog-to-livestock transmission. It is estimated that 7096 of the 
total dog population in Addis Ababa is stray dogs and rabies transmission is more likely in 
stray dogs than owned dogs. In most cases a stray dog which bites a human will run away. 
This complicates the quarantine process and the stray dog continues to spread the virus. 
Fig.7 shows rabies prevalence in the livestock population within and around Addis Ababa. 
According to the information we got from EPHI there is no rabies vaccine for livestock being 
used in Ethiopia and little is known about the status of rabies in livestock populations. In our 
simulation we show how rabies vaccination of livestock has the potential to reduce rabies 
livestock infection, if this intervention was applied in Ethiopia. 

The broken curve in Fig.7 shows the data reported from EPHI for rabies infected livestock 
populations. The smooth curve is the simulation of our system. The number of infected 
livestock decreases from 2008 to 2011 and it increases rapidly from 2011 onwards. The 
transmission of the rabies virus in livestock populations is underestimated; it is likely there 
are more livestock infection cases which are not reported to EPHI. 
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Figure 7: Comparison between the reported data and the model simulation of infected 
livestock population from 2008 to 2014. 
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Figure 8: The trend of livestock populations with different vaccination coverage in 40 years 
time. 

Fig.8 shows that the rabies virus increases rapidly in the livestock population for the coming 
fifteen years and will peak in 2026. Aside from control measures which should be taken in 
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dog populations there is a possibility to reduce the disease by applying 25% vaccination 
coverage in the livestock population, which can cut the number of livestock infections 
predicted around 2026 by half. Our simulation shows there is no way to eliminate rabies 
disease by increasing vaccination of livestock alone. With each increase of 2596 vaccination 
coverage the total impact in terms of reductions of rabies infected livestock gets smaller. 
Even 7596 vaccination coverage does not lead to the elimination of the disease. For this 
reason, we recommend that reducing rabies transmission in the dog population is the best 
method for controlling the transmission in human and livestock populations. 
5. Conclusion 

In this paper we have formulated an SEIR model of rabies with vaccination. We have 
computed the basic reproduction number R, «2 and the effective reproduction number 


R, «1.6, which shows rabies is endemic within and around Addis Ababa. The current 10% 


vaccination coverage for dog population is very less. Immunization of dogs have to be 
increased by combining with other intervention strategies. Combination of interventions in 
dog populations can greatly reduce rabies infection in human and livestock populations. A 60% 
vaccination coverage in the dog population reduces the threshold parameter R, less than 
unity, which means rabies will die out from the community. A more aggressive intervention 
consisting of 6096 vaccination coverage in dog populations, 15% culling of stray dogs and 
reducing the annual crop of newborn puppies by 2596 will reduce the number of human and 
livestock infection by 7096 which leads to a fast eradication of the disease from the 
community. In Fig. 9 the upper curves show human and livestock infection cases before 
interventions in dog populations. After combining interventions in dog populations, the 
number of infected cases in livestock and humans are greatly reduced as shown by lower 
line curves. 

According to the sensitivity analysis that we have performed the annual birth of dog 
populations and dog-to-dog transmission rate are the most positively sensitive parameters. 
The natural death rate of dog and disease induced death rate are the most negatively 
sensitive parameters. Using numerical simulation of our model current data predicts that 
rabies transmission within and around Addis Ababa will increase in the coming years. It will 
reach to the peak in 2024 by killing more than 150 people and finally it will become stable. 
The number of annual birth of newborn puppies greatly affects rabies infection in human 
and livestock. As we decrease the annual birth of puppies by 3096, human and livestock 
infection also decreases by the same percent. Therefore, controlling the annual birth of 
newborn puppies is one of the best ways to reduce human and livestock rabies infection. 
Rabies transmission in livestock population increases rapidly for the coming eighteen years. 
It will reach to the peak in 2024 by killing more than 180 of them. In our simulation we have 
applied vaccination for different number of livestock populations. With 2596 vaccination 
coverage in livestock population, we can save half of them that will be dying from rabies in 
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the coming thirty years. This can be observed from Fig.10 by taking the difference between 
the curves of no vaccination and 25% vaccination. 
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Figure 9: Effects ofcombining interventions in dog populations (cid) on human and livestock 
rabies infection. 

As we increase vaccination coverage for livestock populations the number of livestock 
infection reduces, we strongly recommend that more control measures should be adopted 
to save many numbers of vulnerable livestock populations. But this should be combined with 
the control measures which we recommended for dog population. 
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Figure 10: Comparison between no vaccination and 25% vaccination for rabies infected 
livestock. 
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